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The Eigenvalue Analysis of the Density Matrix of 4D Spin Glasses Supports Replica 

Symmetry Breaking 
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We present a general and powerful numerical method useful to study the density matrix of spin 
models. We apply the method to finite dimensional spin glasses, and we analyze in detail the 
four dimensional Ed wards- Anderson model with Gaussian quenched random couplings. Our results 
clearly support the existence of replica symmetry breaking in the thermodynamical limit. 
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I. INTRODUCTION 

Replica Symmetry Breaking (RSB)@ was introduced 
more than twenty years agoEl as a crucial tool to describe 
the low temperature phase of spin glassesu. One can 
see replicas as an extension of Statistical Mechanics that 
can be very useful wjijen studying complex systems, such 
as structural glassesu or spin glassesB, where the ergod- 
icity breaking in the low temperature phase cannot be 
described with the help of an infinitesimal external con- 
stant magnetic field. 

If on one side there is little doubul left about the cor- 
rectness of the RSB description of the low temperature 
phase of thtt-jaiGan field models, on the other side the 
controversyoQala regarding its applicability to finite di- 
mensional systems such as realistic, physical spin glasses, 
is alive and in good health. 

Unfortunately, we are only starting to guess how to 
address the question of the existence of RSB in reaJ spin 
glasses from a truly experimental point of viewt£l: be- 
cause of that, and because of the inherent very high com- 
plexity of the relevant analytic computations, most of the 
recent progresses are coming from numerical simulations. 

The output data of numerical simulations are never 
as reliable as analytic (and, even better, rigorous) re- 
sults. So if on one side the results of .auincrical simu- 
lations of four dimensional spin glassea§lLil support the 
RSB scenario (as indeed happens for the three dimen- 
sional model!) on the other side one can argue that these 
indications could turn out to be fallacious on larger lat- 
tices, on longer time scales, at lower temperatures... (see 
for example |l2| for a typical criticism to typical numerical 
simulations). 

It is clear that new approaches to this important issue 
are-precious: Sinova, Canright, Castillo and MacDon- 
aldlij have recently proposed such a new tool, that can 
allow a better study of spin glasses. They have noticed 
that the spin-spin correlation matrix {(JiCFj) (that we will 
discuss in better detail in the next section) shares the 
main-, properties of a quantum mechanical density ma- 
trixes: it enjoys positivity, hermiticity and has unit trace 
(notice that our normalization differs from theirs, see 



reversal symmetry is broken, and thus one should expect 
at least one non vanishing eigenvalue of the density ma- 
trix in the thermodynamical limit, due to the extended 
charactcr_of the eigenvector related jsrith the symmetry 
breakingES What is new is the claintJ that the presence 
of RSB is equivalent to the existence of more than one 
non vanishing eigenvalues in the thermodynamic limit. 
Thus armed the authors of [13] undertook the study of 
the Edwards- Anderson model with Gaussian couplings 
in four dimensions, were they found results that they 
judged inconsistent with the detection of RSB on lattices 
of linear size up to 6 (i.e. of volume up to 6 ). 

The efforts of |l3|were limited to such small lattice sizes, 
because the memory and the numerical effort required 
in their approach grows as L 2D (in the following L will 
be the lattice linear dimension, and D the space dimen- 
sionality). It is clear that their simulation strategy and 
data analysis-can sometimes go wrong, as it is evidenced 
by its failureEj in the analysis of the Random Field Ising 
model. In that case, only turning to the standard numer- 
ical strategy^, that focuses on the Raijisi order parameter 
function, P(q), they could establish^ the (plausible) ab- 
sence of RSB in this model. 

Here we present a numerical strategy for the study of 
the density matrix of spin glass with a cost of the or- 
der L D . We propose a more convenient data analysis, 
given the expected behavior of the density of eigenvalues 
of the density matrix in the thermodynamic limit (see 
next section and reference |l6]) . In this way we have been 
able to study the Edwards- Anderson model with Gaus- 
sian couplings on lattices of volume up to 8 4 , at the same 
temperatures .as in [l3| We obtain results that support an 
RSB scenarios. Very interesting information about the 
density matrix in a RSB scenario can also be obtained 
through mean field calculationstll Moreover the numeri- 
cal approach that we have developed here can be applied 
to any spin model. 

When completing this manuscript, a note reporting 
another efficient approach Jjp the density matrix spec- 
tral problem has appearedO. In this work Hukushima 
and Iba deal with the four dimensional spin glass model 
with binary (rather than Gaussian like in our case) cou- 
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10 4 , reaching the same conclusion that we present here, 
i.e. arguing for the presence of RSB in the infinite vol- 
ume limit (they also discuss an interesting method for 
studying temperature chaos). 

The layout of the rest of this paper is as follows. In 
section II we define the model and the associated density 
matrix, discussing its basic properties and the numeri- 
cal approach of |l3|. Our own strategy is presented in 
subsection y A| , and a working example is analyzed in 
subsection pi B| , where the (replica symmetric) ferromag- 
netic Ising model in four dimensions is analyzed. Our 
numerical simulations of the Edwards-Anderson model 
in 4 dimensions are described in section III, Our results 



are presented and discussed in section IV. Finally, we 
present our conclusions in section [v|. 



attention to this quantity. Reference [T| wisely suggested 
to look at the correlation function of a single sample as 



a matrix, Cij. We define here Cij as 



L D 



(5) 



(notice the difference in the factor L~ D with the def- 
inition of references [D|[l5|). The gauge transformation 
(^) acts on the matrix Cjj as an unitary transformation. 
Therefore, contrary to the individual elements of Cjj it- 
self, the spectrum of Cjj does not become trivial after the 
disorder average. It is easy to checka that Cjj is sym- 
metric, positive definite, and has trace equal to one, just 
like a quantum mechanical density matrix. Thus the cor- 
responding eigenvalues, 1 > Ai > A2 > . . . A^ > 0, verify 



II. THE MODEL AND ITS DENSITY MATRIX 

We consider the four dimensional Edwards-Anderson 
spin glass in a periodic box of side L. The N elementary 
spins can take binary values, — ±1, and they are de- 
fined on the vertices of a single hyper-cubic lattice of size 
V = L D . We consider a first neighbor interaction: 



H 



(id) 



(1) 



The quenched couplings, Jij = Jjj, are drawn from a 
symmetric probability distribution function of zero aver- 
age and variance J 2 . It is customary to take J as unit 
of temperature, and then to set J = 1 : this is what we 
do. Two popular choices are the one of a binary proba- 
bility distribution Jjj = ±1 or to take J Gaussian dis- 
tributed. Here, we draw the quenched random couplings 
from a Gaussian distribution (also in order to allow a di- 
rect comparison with the work of ^3j) . For all the relevant 
observables one first compute the thermal average on a 
single realization of the couplings (sample), hereafter de- 
noted by (...), and later the average with respect to the 
couplings is performed (we denote this disorder average 
by an ovcp,linc). The model ([j]) undergoes a spin glass 
transition^ at T c = 1.80 ± 0.01. 

The average os^er the couplings Jij induces a (trivial) 
gauge invarianceO in the model. If one chooses a generic 
binary value for each lattice site, r\i = ±1, disorder aver- 
aged quantities are invariant under the transformation 



J, 



0", 



(2) 
(3) 



Now let iji be a random number that takes the probability 
i the values ±1. If one considers the spin-spin correla- 
tion function, the symmetry (||) yields the disappointing 
result that 



(piOj) = rjiijjiaicrj) = S itj , 



(4) 



(that is true since this relation is valid for every value of 



fc=i 



(6) 



Following [TJ the authors of |l3| have argued that in the 
paramagnetic phase all the A& are of order j^, and thus 
vanish in the thermodynamical limit. On the other hand 
in the spin glass phase time reversal symmetry is bro- 
ken, which implies some non local ordering pattern for 
the spins (unfortunately only known by the spins them- 
selves), and hence at least one eigenvalue, Ai should re- 
main of order one when N — > 00. They also claimed that 
presence of RSB is equivalent to more than one eigenvalue 
being of order O(N ) when N — > 00. Furthermore they 
stated that each non vanishing eigenvalue corresponds to 
a pair of pure states: the correspondence to a pair of pure 
states is because of the global a — ► — a symmetry of the 
Hamiltonian ([j]) and of the matrix Cjj. Notice that this 
might be a clue for the solution of the formidable pral 
lem of defining pure states in a finite volume systemLH 
The fact that the presence of more than one extensive 
eigenvalue (of order O(N )) when N — > 00 is equivalent 
to RSB is true in the mean field picture, as can be veri- 
fied in a_mean field analytic computation at the first step 
of RSBE3. 

Combining perturbation theory and droplets ideas it 
was also possible to tellL2l that in a non RSB scenario the 
second eigenvalue should not decay slower than 



A 2 ~ L~ 



where the droplet exponent in four dimensions is 
0.6-0.8. Actually when the lattice size is larger than the 
correlation length (which might Jiot be the case in the 
achievable numerical simulationgI3) , they expect a much 
faster decay. 

Using_the parallel tempering optimized Monte Carlo 
schemeO, the authors of [II] calculated the matrix Ci.j, 
(a computational task of the order L 2D , since the lack 
of translational invariance prevents the use of the Fast 
Fourier transform). They eventually diagonalized the 
matrix. When comparing results for different disorder 
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Afe, that they tried to characterize by their mean and 
typical value. They found that the mean and the typi- 
cal value of the second eigenvalue were decreasing as a 
function of lattice size in a double logarithmic plot for 
lattices up to 6 4 (see figure 7 of the second of references 
|l3|) . Because of that they argued about the absence of 
RSB in the model. 



A. An Effective Approach to the Study of the 
Density Matrix 

Studying the spin-spin correlation function aj by an- 
alyzing the usual density of states g u 



replicas a; and a 7 : 



1 N 



(8) 



k=l 



would not work: because of the constraint (J6|) in the 
TV — > oo limit g u (X) is a normalized distribution function 
with support in the [0, 1] interval with mean value 0. 
In other words, this definition implies that in presence 
of a generic finite number of extensive eigenvalues for 
large volumes g u (X) = <5(A), which does not contain much 
information. 

In our case we cannot weight all the eigenvalues with 
the same weight: to consider a sensible indicator we can 
decide to use as weight A& itself, and to define the mod- 
ified density of states of the matrix c,^: 



(9) 



k=l 



It is natural to expect g(X) to converge in the N —> oo 
limit to a function containing a continuous part, plus a 
delta function at A = (because a number of order O(N) 
of eigenvalues will-be of order 0(N~ 1 )). A calculation 
at one step of RSBt£l tells us that this is indeed the case. 
Moreover in the one step calculation the continuous part 
do not show any gap, and it covers all the interval be- 
tween A = and A = <7ea, the Edwards- Anderson order 
parameter (see also figure 1 of the second of reference 13). 
Therefore, from the point of view of checking Replica 
Symmetry Breaking, to concentrate on the behavior of 
individual eigenvalues does not look the best strategy. 
Instead, as it is customary when analyzing density of 
statesEa, one can start by considering the moments for a 
single disorder realization, gj{X): 



dXX r gj(X) = 



N 

=Trc r+1 

fc=i 



(10) 



Our main observation is that we can compute the trace of 
the r-th power of the matrix c, using r real replicas (inde- 
pendent systems, with the same realizations of quenched 



i=i 



Then it is easy to show that 
Trc r = (q Ql < Q V 



■Q 



(11) 



(12) 



Thus, for instance, the (disconnected) spin glass suscepti- 
bility is xsg = N Tr c 2 . In this language the relationship 
between non vanishing eigenvalues and the phase tran- 
sition from the paramagnetic to the spin glass phase is 
very direct. 

It is now very easy to suggest a numerical strategy of 
order L D : make the Monte Carlo simulation in parallel 
for a discrete number of replicas, and use them to calcu- 
late the appropriate number of moments of gj(X). Then 
use this information to extract the largest eigenvalues of 
the matrix c. Unfortunately standard methods foe ex- 
tracting the probability density from its momentsE3 use 
orthogonal polynomials. Clearly, given the limited nu- 
merical accuracy that we can expect to obtain for the 
Tr c r , the use of orthogonality methods is out of the ques- 
tion. We have instead used a cruder method. We define 
a cost function 



,60 



E 



EI 



=i £fc 



Trc' 



(13) 



and minimize it, using the values of the at the min- 
imum as an approximation to the eigenvalues. This 
method can be checked on small lattices, using the di- 
rect computati on of c and its eigenvalues. It turns out 
(see subsection [I B and section |EV| ) that it is extremely 
precise for the first eigenvalue, Ai, but that already for 
the second eigenvalue, A2 the systematic error is at the 
10% level using 12 replicas. Fortunately we can do better 
than setting A2 ~ £2- Let us define a (further) modified 
density of states in which we do not include the first 
eigenvalue 



N 



g(X) = A fc S(X - A fc ) . 



fc=2 



Its moments are 



11 



dXX r g(X) = [Trc r +! 



(14) 



(15) 



where we have denoted by A r the subtracted traces. The 
right hand side of equation (|l5|) can be accurately calcu- 
lated using the cost function, and contains all the infor- 
mation that we need. 

One could still worry about the bias induced by our 
use of the cost function to obtain Ai. This can be easily 
controlled, because, since the eigenvalues of the matrix 
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a situation where we can expect that A r is clearly and 
substantially larger than A r+ i. On the other hand, if the 
bias on Ai is 8, it will affect A r of a quantity of the order 
of (S r Ap 1 ). Therefore, a bias dominated subtracted 
trace will be characterized by succ essiv e moments of g(X) 
being very similar (see subsection [IB). 

Let us conclude this subsection by discussing the dif- 
ferent scenarios that could describe the scaling of the 
subtracted traces, in the L — » oo limit. For a standard 
replica symmetric model, such as the usual ferromagnetic 
Ising model, we expect A r+1 = 0{L ~ rD ) . In a RSB 
scenario we expect that for L — > oo A r+ i tends to a 
finite value (and that finite volume corrections due to 
the eigenvalues that create the S(X) in g(X) are of the 
form 0(L~ rD ), while iother finite size corrections due 
to critical fluctuations may not decay so fast). Finally, 
in a droplet scenario, if one assumes that the subtracted 
traces are controlled by A2 , then equation (^) implies that 

A7 = O(L~ r0 ) , (16) 

with 9 = 0.6-0.8 in four dimensions (recall that this is 
an upper bound in the decay of A2). The only way out 
from this scaling behavior in a droplet picture would be 
to assume that a number of the order ifi (£ > 0) of 
eigenvalues, is of order L~ B : we are not aware of any 
argumentii3 that would imply the existence of a divergent 
number of critical eigenvalues in a droplet picture. 



B. A Simple Example: the Ferromagnetic Ising 
Model 

As a first check we have studied the ferromagnetic Ising 
model in four dimensions. Here the Hamiltonian has the 
same form than in (|l|), but with Jy = 1. We have stud- 
ied the system at T = 0.5T C , to prove the deep broken 
phase with small correlation length (the critical temper- 
ature is hereB T c = 6.68025 ± 0.00004)). We have simu- 
lated in parallel (in this case without parallel tempering, 
but with an usual heat-bath updating scheme) twelve 
replicas of lattices of linear size L = 3,4,6 and 8, for 
3 x 10 5 Monte Carlo steps, starting from a fully ordered 
state. 

In this simple case the density matrix cy can be very 
easily diagonalized. The correlation function {oi<Jj) de- 
pends only on the distance between the two spins, Xi — Xj, 
and thus the eigenvectors are proportional to exp[ifc • a?,-], 
where the wave-vectors k verify the usual quantization 
rules on a periodic box. It is straightforward to show 
that the corresponding eigenvalues are 



h = 



N 

E 



L D 



(17) 



and, given the ferromagnetic character of the interaction, 
the largest eigenvalue corresponds to k = (the magne- 
tization, M): 
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FIG. 1: Cost function ( |l3| ) estimate of the largest eigenvalue 
of the density matrix, as a function of the number of calcu- 
lated moments (see equation (E^)), for the four dimensional 
Ising model at T = 0.5T C , in a L — 4 lattice. The horizontal 
lines correspond to {M 2 } plus or minus one standard devia- 
tion. 



0.96114 



0.96112 
0.96110 
0.96108 
0.96106 
0.96104 
0.96102 
0.96100 



ISING FM 



L = 8 
m 2 = 0.961133 



6 7 8 9 
# of moments 



10 11 



FIG. 2: As in figure [j] but for a L = 8 lattice. 

In figures |l| and || we compare our estimate of Ai for 
the L = 4 and L — 8 lattices, as obtained from the 
magnetization (the horizontal band is (M 2 ) plus/minus 
an standard deviation), and from the cost function (13|). 
As both figures show, 12 replicas are surely enough to 
obtain agreement within errors, which in this case are 
particularly small. 

Having gained confidence in our procedure we can now 
check evolution of the subtracted traces with increasing 
lattice size (figures || and ||). The two values are very 
small, decreasing with the lattice size and almost (but 
not completely) compatible with zero. One should no- 
tice that A3 and A2 are compatibles within errors for all 
lattice sizes (we will see in section |^ that in the spin 
glass case the situation is very different): in the ferro- 
magnetic case the real A 3 and A 2 are so small that they 
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FIG. 3: The subtracted trace, A2, as a function of the lattice 
size, for the four dimensional ferromagnetic Ising model. 
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60% level. For each replica we have measured the perma- 
nence histogram at each temperature, and we checked its 
flatness. We controlled thermalization by checking that 
there was no residual temporal evolution in the Binder 
cumulant and in Trc . 

The main scope of the simulation has been to obtain 
Trc'', for r = 2, . . . , 12, using equation There is 

an awfully large number of equivalent ways of forming 
the trace q ai ' a2 q a2,a3 . . . q a r,<n w hen one may choose the 
replica labels dj out of twelve possible values. One needs 
to find a compromise between loosing statistics and wast- 
ing too much time in a given disorder realization (the dis- 
order average is the critical factor controlling statistical 
error). Our compromise has been the following: given 
the special importance of this observablefl we have calcu- 
lated the 12 ( 1 ^~ 1 ) possible overlaps g Ql >° 2 , and we have 
computed Tr(c 2 ) using all the 66 quantities. For traces of 
higher order we have considered only twelve contributions 
of the form g M+i g i+i,i+2 _ _ _ q i+r,i^ for i = 1, 2, . . . , 12 
(the sums are understood modulo 12). 
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FIG. 4: As in figure a but for A 3 
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previous subsection. One might ask how come that we 
were able to resolve such an small bias, given the com- 
paratively large errors reported in figures [j] and g: this is 
due to the strong statistical correlations betweenTr(c r ) 
and our estimate for X[. 



III. THE MONTE CARLO SIMULATION 

We have studied by numerical simulations the four di- 
mensional Edwards- Anderson spin glass with quenched 
random Gaussian couplings ([j]). We have simulated 12 
real replicas in parallel using a heath bath algorithm and 
Parallel Tempering^, on lattices of volume 3 4 , 4 4 , 6 4 and 
8 4 . The ratio between full lattice heat bath sweeps and 
parallel tempering temperature swap attempt was one 
to one. For all lattice sizes the largest temperature was 
Tmax = 2.7 and the lowest temperature T m ; n = 0.8 (see 
table I for details of the numerical simulation). The prob- 

-l _r ±1 _ i i l ± _ j_ _ 



FIG. 5: The Binder cumulant as a function of temperature, 
for the 4D Edwards- Anderson model on lattices of linear size 
L = 3, 4, 6 and 8. 

In addition to the Tr(c r ) we have measured the Binder 
cumulant (see figure ||) We have also measured a second 
adimcnsional operator 



Trc 25 

that we show in figure ^. .— . 

The theory of finite size scalingcil predicts that adi- 
mensional quantities close to criticality are functions of 
where v is the thermal critical exponent (in 
D = 4 one findsE3J v = 1.0 ± 0.01)). The crossing points 
signals the spin glass transition at T c — 1.8 with similar 
accuracy for both the cumulants that we have consid- 
ered. At the lowest temperature that we have reached 
the L = 6 and L = 8 lattices seem to be far enough from 
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FIG. 6: As in figure ||, but for the Bi cumulant. 



IV. NUMERICAL RESULTS 

To compare our results with the ones of [l3] we will spe- 
cialize here to T — 1.0. We start by checking on small 
lattice sizes (see in figure the L = 4 data) the cost func- 
tion procedure. In this case the estimate of Ai that one 
can obtain by using the cost function can be compared 
directly with the result obtained by diagonalization of c: 
we find a fair agreement. For larger lattices we can only 
check the convergence of A,, as a function of the number 
of moments (see figure ||). Again, the convergence looks 
fast enough for our purposes. We show in figure ^ the 
probability distribution of Ai . The low eigenvalues tail 
is basically lattice size independent. 

We show our results for A2 and A3 in figure [l(] and 
figure [n], respectively. A2 is a factor of 10 larger than 
A3: our dat a are not bias dominated (see subsections 
II A| and |II B| ) . The fact that the data point for A3 in the 
L = 8 lattice is above the L = 6 one and at two standard 
fluctuations from compatibility may be due either to a 
strong fluctuation, or to a first glimpse of bias effects. If 
one sticks to the bias hypothesis, the effect on A 2 can be 
(very conservatively) estimated as the difference of the 
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FIG. 7: Disorder averaged cost function ([[3]) estimate of the 
largest eigenvalue of the density matrix, as a function of the 
number of calculated moments (see equation ([TJ])), for the 
four dimensional Edwards- Anderson spin glass at T = 1.0, 
on a L = 4 lattice. The horizontal lines correspond to a 
numerical diagonalization of the matrix Cjj plus or minus a 
standard deviation. 
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FIG. 8: Disorder averaged A2 for the four dimensional 
Ed wards- Anderson spin glass at T = 1.0 as a function of 
the number of computed moments, on different lattice sizes. 
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■^V samples 


N measures 


-^thermal 


Np 


3 


2800 


50000 


50000 


20 


4 


2800 


50000 


50000 


20 


6 


1208 


150000 


150000 


40 


8 


362 


100000 


200000 


40 



TABLE I: Relevant parameters of the Monte Carlo simula- 
tion. L is the lattice size. A r samp i ea denotes the number of 
realizations of the Gaussian couplings. The number of Monte 
Carlo steps (heat bath sweep plus temperature swap attempt) 
discarded for thermalization was iVthermai- Np is the number 
of temperatures simulated in the parallel tempering. Finally, 
measures were taken during' AT„„, Monte Carlo stens. 



L = 6 and L = 8 data points corresponding to A3. This 
difference is well covered by the error in the L = 8 data 
point for A2. 

After the above considerations we can now proceed to 
the infinite volume extrapolation. In figure [l^ we plot 
the data for A2 as a function of L~ D . It is evident that, 
letting aside the L = 3 data, a linear fit is appropriate. 
The extrapolation to infinite L is definitely different from 
zero: 



L > 3, A 2 = 0.0119 ±0.0003 
L > 4, A2" = 0.0102 ± 0.0004 



X 2 /dof = 17.8 , (20) 
X 2 /dof = 1.73. (21) 



In figure we plot the data as they should scale accord- 
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FIG. 9: Probability distribution of the largest eigenvalue as 
calculated in the four dimensional Edwards- Anderson spin 
glass at T = 1.0, for lattices of linear size L — 4, 6 and 8. 
The binning in the L = 8 lattice was reduced by a factor of 
two, due to the smaller number of samples. 
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FIG. 10: Disorder averaged subtracted trace A2 for the 
four dimensional Edwards- Anderson spin glass at tempera- 
ture T = 1.0 as a function of the lattice size. 



equation (|T^) yields a very high value of % 2 /dof either 
when we include the L = 3 data or when we-exclude 
them (we use 9 = 0.6, the lowest possible valueEJ): 



L > 3 , x7 dof 
L > 4 , x 2 /dof 



17, 
14. 



(22) 
(23) 



CONCLUSIONS 



We have proposed and used a new numerical approach 
to the study of the density matrix in spin glasses. The 
original idea of namely to introduce the density ma- 
trix in the spin glasses context, allows to make interesting 
calculational3, and might even proyejiseful to the defini- 
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FIG. 12: Disorder averaged A2 as a function of L~ for the 
four dimensional Edwards- Anderson spin glass at T = 1.0. 
The dashed line is for a linear best fit, excluding the L — 3 
data. 



Our method is a further step beyond the useful ap- 
proach of of [l3|. The technology we have developed can 
be safely applied to the study of different spin models. 
The main limitation of our approach is not related with 
the use of the density matrix, but with the extreme diffi- 
culty in thermalizing large lattices deep in the spin glass 
phase. Should an efficient Monte Carlo algorithm be dis- 
covered, our method would be immediately available, be- 
cause the computational burden grows only as L D . Very 
recently, another optimized method has been proposed 
by Hukushima and IbaEZl. Using their method they were 
able to study 10 4 lattices, using binary rather than Gaus- 
sian couplings (which strongly speeds up the simulation). 

Using our approach we have been able to show that 
the density matrix approach for the four dimensional Ed- 
wards Anderson model with Gaussian couplings in lat- 
tices up to L = 8, and temperatures down to T = 1.0 

/ n rnm \ • r__^^__ • . . -ii t-» ciTi •_ j J 
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predicted by the alternative droplet model. In this re- 
spect, the results are in full agreement with the avail- 
ables studiesa of the Parisi order parameter, and with 
the recent results of [l7| A word of caution is in order: 
the (postulated) impossibility of getting thermodynamic 
data in the reachable lattices sizeslL3, affects equally to 
the P(q) approach and to the density matrix approach. 
However our data for adimensional quantities, such as 
the Binder or B3 cumulant, seem very hard to reconcile 
with the possibility of a purely finite volume effect. 
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FIG. 13: Disorder averaged A2, as a function of L~ for the 
four dimensional Edwards- Anderson spin glass at T — 1.0. 
The droplet 9 exponent is chosen at its lower bound, 9 = 0.6. 
The dashed (dotted) line is for a linear best fit, excluding 
(including) the L = 3 data point. 



that there are serious difficulties with the scaling laws 
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